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Column abundance estimates of sodium, and analogously, 
potassium, in Mercury’s exosphere are strongly correlated to the 
surface reflection model used to calibrate the spectral data and 
the surface reflection model incorporated into the atmospheric 
radiative transfer solution. Depending on the surface reflection 
model parameters used, there can be differences in calibration 
factors of up to ±30% and differences in estimated column abun- 
dance of up to ±35%. Although the surface reflectance may not 
be used in the calibration of spacecraft measurements, the inter- 
action between the reflected surface light and the atmospheric 
brightness remains important. < 1997 Academic press 


INTRODUCTION 

The sources and sinks of Mercury’s atmospheric sodium 
and potassium are not well known. Efforts to understand 
them better require knowledge of the distribution of Na 
and K about the planet. Linking spatial variations in col- 
umn abundances to potential sources and sinks would in- 
crease our knowledge of Mercury’s surface and atmo- 
spheric chemistry. This study examines the effect of chosen 
parameter values for the surface reflectance model on the 
calculation of mercurian atmospheric column abundances. 
We show that the reflectance properties of the mercurian 
surface can have significant effects on the abundances cal- 
culated from spectroscopic data. The results of this study 
define what level of column abundance spatial variability 
could be ascribed to current data reduction methods. 

We assume that groundbased observations of spectral 
lines are calibrated with reference to the brightness of 
Mercury’s surface at a nearby wavelength. The advantage 
of this procedure is that the calibration data are from the 
same exposure as the line data. They share all the sky 


conditions, such as clouds near the horizon. The disadvan- 
tage, one of the points emphasized in this paper, is that 
the brightness of the surface is somewhat uncertain. Cali- 
bration against a standard star may be preferable for obser- 
vations from a spacecraft, but not from the ground. 

Regardless of the methodology used for the calibration 
and column abundance determinations, the surface reflec- 
tion theory used in the mercurian atmosphere literature is 
Hapke’s model (Hapke 1981, 1984. 1986). Our study focuses 
on the effects of the choice of Hapke model parameters on 
the calculated column abundances of Na in Mercury's atmo- 
sphere. Here we work with six Hapke parameter sets, each 
of which satisfies the disk integrated photometry of Danjon 
(1949) for Mercury. To accomplish our goal we have re- 
duced and analyzed sodium emission measurements from 
J anuary 23, 1 988 (Sprague etal 1 996a), obtained at the Cata- 
lina Observatory outside Tucson, Arizona. The observa- 
tions were made using the LPL echelle spectrograph, 
mounted in the focal plane of the 1.5-m telescope. 

ROLE OF HAPKE’S MODEL 

Hapke’s model (Hapke 1981, 1984, 1986) is a function 
of incident, emission, and phase angles. The parameters 
include the single scattering albedo, n\ the opposition ef- 
fect parameters b 0 (amplitude) and h (width), the surface 
roughness parameter, 0 , and the single particle scattering 
function (which may have one to three parameters). Table 
I and Figure 1 display six viable solution sets to the 
Hapke model parameters based on modeling Danjon's 
(1949) disk-integrated phase curve of Mercury. The first 
three Hapke model solutions listed use a two parameter 
Legendre polynomial of the form 

p(a) = 1 + b cos(a) T c(3 cos 2 (a) - 1 )/2 
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TABLE I 

Hapke Parameter Solution Sets to Mercury’s Disk-Integrated 
Phase Curve 


Parameter 

Set 1" 

Set 2“ 

Set 3" 

Set 4 

Set 5 

Set 6 

w 

0.21 

0.20 

0.23 

0.22 

0.218 

0.286 

bo 

1.85 

2.4 

2.5 

1.0 

1.0 

1.0 

h 

0.03 

0.11 

0.09 

0.0604 

0.0284 

0.0438 

b 

0.40 

0.20 

0.18 

0.268 

0.266 

0.307 

c 

0.40 

0.18 

0.15 

0.276 

0.148 

0.531 

0-bar 

20 

21 

25 

22 

20 

32 

rms 

0.121 

0.162 

0.212 

0.118 

0.124 

0.117 


" Used a Legendre polynomial for the single particle scattering function. 
Others used a double Henyey-Greenstein function. Error bars: w = 
±0.02, bo = ±0.02, h = ±0.01, b = ±0.05, c — ±0.05, 0-bar -- ±5. rms 
is the root mean square difference between the data and the model. It 
is defined as the square root of the sum of: the differences squared divided 
by the number of data points. 


for the single scattering function, where a is the phase 
angle. The first solution is taken from Bowell et al (1989) 
and the second and third solutions are solutions 1 and 2 


from Veverka et al (1988). The last three Hapke model 
solutions use a two parameter double Henyey-Greenstein 
function of the form 

p(a)= + <0 - b 3 

(1 - 2b cos(a) + b 2 y 2 (1 + lb cos(a) 4- b 2 ) 312 

for the single scattering function, where a is the phase 
angle. These solutions were found using a modified least 
squares grid search and are published here for the first time. 

Surface reflection modeling enters into the calculation 
of sodium and potassium column abundances in two ways: 
(1) in determining the instrument sensitivity and (2) in 
calculating the emission from atoms induced by the ab- 
sorbtion of reflected light from the surface. Absolute cali- 
bration assigns a value of brightness in rayleighs to the data 
numbers measured by the imaging device. The method of 
absolute calibration for atmospheric emissions from Mer- 
cury’s atmosphere is reviewed in the next section; however, 
the details are available in Sprague et al. (1996b). 

Calculations of column abundances must include reso- 
nant scattering induced by surface reflected light. Two dif- 



FIG. 1. Disk-integrated solar phase curve for Mercury. The open diamonds are Danjon’s (1949) V filter observations of Mercury. The values 
of the parameter solution sets are listed in Table I. 
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TABLE II 

Observing Parameters 


Date UT 

23 January 1988 

Frame No. 

168 

UT start time 

21 hr 10 min 

Exposure time 

50 sec 

Mercury-Sun distance 

0.3239273 AU 

Phase angle 

70.8° 

Diameter 

6.45 arcsec 

Sub-Earth longitude 

70.6° 

Subsolar longitude 

0.0° 

Mercury-Sun doppler shift 

-141 mA 

Solar flux at Mercury 

4.965 X 10 14 quanta/cm 2 sec A 

g-factors (D2, Dl) 

36, 20 sec 1 


ferent approaches have been used to calculate the radiative 
transfer of sodium emission at Mercury. Killen (1988) and 
Killen et al (1990) have used Chandrasekhar’s solution to 
the planetary problem in terms of X and Y functions and 
their moments. The monochromatic results are then inte- 
grated over the lineshape. Hunten and Wallace (1993) 
adopted a doubling-adding numerical method with the ap- 
proximation of complete frequency redistribution (CFR). 
Though the results of the two methods have not been 
compared, both should have ample accuracy. Both models 
assume a Maxwell-Boltzmann distribution for the velocity, 
which may not be the case for Mercury's atmosphere. How- 
ever, for this study we are concerned only with uncertain- 
ties related to the surface reflection model. 

ALGORITHM FOR TESTING DATA 
REDUCTION METHOD 

Using an actual spectrum showing Na emission along a 
north-south (N/S) slit across Mercury’s Earth-facing disk 
(Sprague et al. 1996a), we calculated the column abun- 
dances six times, once for every set of Hapke model param- 
eter solutions given in the previous section. The goal was to 
examine any changes in the estimated column abundances 
owing to the choice of Hapke parameters. 

The Atmospheric Sodium Emission Data 

Table II gives the physical parameters of Mercury for 
the atmospheric measurements considered here. For this 
particular observation the spectrograph slit was placed 
across the illuminated sector, in the N/S direction, of a 
gibbous planet (70.8° phase angle) and offset from the 
central meridian by —2 arcsec or —0.6 mercurian radii. 
The N/S data scan is divided into five sectors which are 
correlated with a physical position (latitude and longitude) 
by using an appropriate model which includes the effects 
of seeing smear (<x). The seeing a is defined as the half 
width at the Me point of a two-dimensional gaussian. Esti- 


mations of the seeing smear were made by generating a 
model of the surface reflectance and smearing it to a variety 
of seeing a. The seeing <j for the data is determined by 
comparing a N/S scan from each version of the smeared 
model to the actual planetary data within the continuum 
portion of the spectrum and finding the scan from the 
smeared models which most closely approximates the 
shape and relative size dimensions of the continuum scan 
from the data. Table III lists the a values calculated using 
each of the Hapke model parameter solution sets to gener- 
ate the initial surface reflectance model. 

The instrument sensitivity or calibration factor (CF = 
brightness giving 1 DN/sec) is given by 

CF = B S IQ krayleigh sec DN _1 , 

where B$ is the surface brightness and Q is defined as the 
average number of data counts (DN) per angstrom per 
second. The surface brightness (B s ) is given by 

B$ = E\ X 47 tR r - 10" 12 Mrayleigh Al 

where is the solar flux at Mercury at wavelength A and 
R r is the average surface reflectance from the center of a 
swath across the smeared model corresponding to the slit 
across Mercury from which the data were taken. The values 
for 5 S , Q , and CF are listed in Table III as derived from 
each of the six Hapke model parameter solution sets. 

Column Abundance Calculations 

Column abundances were calculated using the method 
of Hunten and Wallace (1993) which accounts for the self- 
shadowing of sodium radiation by its own column of atoms 
above the surface. Calculations of the sodium column 
abundances were made using this method for each of the 
Hapke model parameter solution sets with the correspond- 
ing seeing a and calibration factors. The contribution of 
sodium emission due to the resonant scattering induced 
by the absorbtion of surface reflected light depends directly 


TABLE III 

Data Reduction Parameters 


Solution set Nos. 

(T 

B s 

Q 

CF 

1 

1.74 

48.6 

279.09 

174.1 

2 

1.79 

54.3 

279.09 

194.5 

3 

1.81 

59.4 

279.09 

212.8 

4 

1.79 

54.7 

279.09 

196.2 

5 

1.74 

54.7 

279.09 

196.2 

6 

2.00 

31.8 

279.09 

113.9 


Note. B s is in Mrayleighs A Q is in DN A 'sec CF is in krayleigh 
sec DN 
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upon the calculated surface reflectance values (/? R ) based 
on Hapke's model. In order to separate the Hapke model's 
influence on the calibration factor and seeing a from its 
influence on the atmospheric scattering functions and the 
interpolation and self-shadowing process, we did two dif- 
ferent numerical experiments. First we generated six self- 
consistent models in which the seeing as, calibration fac- 
tors, and suites of column abundances are calculated from 
the same Hapke model parameter values. Second, we used 
the same calibration factor and seeing a (the average of 
the values listed in Table III) for each of six trial runs that 
utilized only the unique atmospheric scattering function 
corresponding to the specific Hapke model parameter set 
in the interpolation. 

RESULTS 

Table III and Figs. 2 and 3 summarize our results. The 
self-consistent trials produced different column abun- 
dances as shown in Fig. 2. Differences of up to ±7% are 
found for the seeing a calculations and differences of up 
to ±27% are found for the calibration factors. The end 
result is differences of up to ±35% in the estimated column 
abundances. Figure 3 shows the results of the second nu- 
merical experiment, where the same values of the seeing 
(rand the calibration factors (1.8 arcsec and 181 krayleigh 
sec DN \ respectively) were used. The results show that 
the choice of parameters used in the Hapke model affect 
the radiative transfer solution up to ±9%, considerably 
less than in the first case. The greatest effect of the Hapke 
model is in the data calibration. 

All of these calculations assume an optically thick atmo- 
sphere. An analogy can be made for an optically thin atmo- 
sphere (as is the case for potassium). The calculations for 
the calibration factor and seeing smear are the same for 
both types of atmosphere; therefore, variations in their 
values due to Hapke parameter values are the same. How- 
ever, the self-shadowing effects are much different. In the 
optically thin case the self-shadowing effects are linear; 
therefore, variations in column abundance estimates due 
to absorbtion by the atmosphere of surface reflected light 
are directly proportional to the variations in the modeled 
surface reflectance values (Rr). In the optically thick case 
the self-shadowing effects are highly non-linear, and there- 
fore have been calculated. 

DISCUSSION 

Calculations of the magnitude of the effect of different 
choices of Hapke model parameters on the estimated prop- 
erties of Mercury's exosphere show that the larger effect 
is in the deduced calibration factor. This will change the 
estimated intensities for the w hole observation. A second, 
smaller, effect occurs due to different estimates of the 


amount of illumination of the atmosphere by surface re- 
flected radiation. 

Hapke model solutions to a disk-integrated phase curve 
are descriptive of the reflectance properties of the average 
surface. If the assumption is made that the surface of Mer- 
cury is homogeneous, then the solution to the disk-inte- 
grated phase curve should provide a reasonably accurate 
description of any subarea on Mercury. Note, however, 
that all six parameter solution sets in Table I are from 
viable fits of Hapke's model to Mercury's disk-integrated 
phase curve. Depending on which of these parameter sets is 
incorporated into the surface reflection model, the column 
abundances may be over- or underestimated by up to 
±35%. In addition, while any variations seen in the column 
abundances as a function of position on the disk are real, 
the magnitude of these variations will also be strongly 
dependent on the choice of parameters used in the surface 
reflection model. If, however, Mercury’s surface is hetero- 
geneous, then the surface reflectance for different areas 
on Mercury are described by different Hapke model pa- 
rameter sets. If this assumption is valid, the calibration 
factor derived for a data set will correspond to a mean 
reflectance for the entire area near the center of the image. 
This reflectance will probably be reasonably well described 
by one of the disk-integrated parameter sets, and the error 
in the calibration factor will be as previously discussed. 
The illumination of the atmosphere by light from the sur- 
face will, however, not be correctly described, and errors 
of up to 5% might be made in the deduced abundance as 
illustrated in Fig. 3. These 5% spatial uncertainties may be 
a lower limit. Murray et al. ( 1 974) found that the brightness 
variations across Mercury’s surface from Mariner 10 data 
are similar to those seen for the Moon, and they estimated 
normal albedos (measured at 5° phase) to range from 0.09 
to 0.45 (Veverka et al. 1988). While 5% uncertainties are 
probably well within the noise of groundbased observa- 
tions, they may be significant for future spacecraft observa- 
tions. 

If Mercury's surface is heterogeneous, then the only 
valid method for determining the surface reflectance prop- 
erties is to model disk-resolved data of the surface. Disk- 
resolved data do not currently exist at the number of view- 
ing geometries needed to adequately constrain the values 
of the Hapke model parameters. Because of seeing smear, 
light along the observing slit in the continuum portion of 
the groundbased spectra does not display spatial variations 
in brightness. Thus the assumption of a homogeneous sur- 
face on Mercury is a valid first approximation considering 
the spatial resolution of the atmospheric data. This will 
not be true for future spacecraft observations, as demon- 
strated by the normal albedo variations seen by Mariner 
10 (Murray et al 1974, Veverka et al. 1988). 

With respect to groundbased observations, we ask which, 
if any, of the Hapke model parameter solution sets pre- 



REFLECTANCE ESTIMATES AND COLUMN ABUNDANCES 


79 


Column Abundance Based on the Sodium D1 Line 




Column Abundance Based on the Sodium D2 Line 
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FIG. 2. The sodium column abundances along a north-south slit of Mercury corresponding to an observation along the illuminated hemisphere 
of a 70.2° solar phase observation. The values of the parameter solution sets are listed in Table I. 


sented here is the most reasonable choice for describing 
Mercury’s surface. This may greatly influence the column 
abundance calculated from the spectroscopic observations. 
One method of determining if one set of disk-integrated 
parameters is a better description of the surface than an- 
other is to compare the fits of each of the parameter solu- 
tion sets to disk-resolved data. Figure 4 shows two scans 
across the disk of Mercury taken from an orange filter 
Mariner 10 image digitized from similar plots shown by 
Bowell et al (1989). The Mariner 10 orange filter has an 
effective wavelength of 575 nm (Hapke et al 1975) com- 
pared to the 550-nm effective wavelength of Danjon’s 
( 1 949) V filter data. Accounting for the differences in wave- 
length and the fact that the spacecraft disk-resolved data 
and the groundbased disk-integrated data were taken by 
two different instruments under extremely different ob- 
serving conditions. Fig. 4 shows that all of the parameter 


solution sets listed in Table I, with the exceptions of sets 
1 and 6, are reasonable fits to the Mariner 10 disk-resolved 
data. The best fits are given by parameter sets 2 and 3 
(also known as solutions 1 and 2 from Veverka et al. 
(1988)). Note that the column abundance estimates be- 
tween the Hapke parameter solution sets 2 and 3 vary by 
approximately ±10%. 

Another consideration is the relationship between the 
parameter values and the physical attributes they describe. 
For example, Hapke’s model assumes that the opposition 
effect is due entirely to shadow-hiding. Therefore b 0 is 
defined as the fraction of wp{a) that is scattered from the 
portion of the particle surface facing the source (Domingue 
and Hapke, 1989). By its physical definition the opposition 
amplitude parameter, should not exceed unity. How- 
ever, Hapke (1990) and Mischenko (1992) have recently 
shown that an additional method for producing the opposi- 
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FIG. 3. The sodium column abundances along a north-south slit of Mercury corresponding to an observation along the illuminated hemisphere 
of a 70.2° solar phase observation. The values of the parameter solution sets are listed in Table 1. Unlike the column abundances plotted in Fig. 2, 
these abundances assume the same seeing tr and calibration factor for all estimates of the column abundance. 


tion effect is coherent backscatter. Hapke et al (1996) 
have shown, through reflectance and polarization studies 
of several lunar soil samples, that the opposition effect on 
the Moon is caused by a combination of shadow-hiding 
and coherent backscatter. Coherent backscatter may also 
produce a portion of Mercury's opposition effect. Values 
of b ty greater than unity may result from this dual opposi- 
tion effect mechanism. 

SUMMARY 

The Hapke parameter values chosen to model the sur- 
face roughness of Mercury (or by analogy, the Moon) affect 
the calibration factor, the estimate of the seeing a, and 
the calculated column abundance. If a large data set is 
analyzed entirely with one parameter set under the as- 


sumption that the surface is relatively homogeneous, then 
relative abundances will not be affected. However, for 
abundances calculated over heterogeneous terrain that is 
not well-described by a single parameter set, relative abun- 
dances will be either under or overestimated. Some condi- 
tions are not suited to the exclusive use of one parameter 
set. For example, if a large, smooth, dark basin is embedded 
in a region with large-scale roughness and high albedo, 
the same Hapke parameter set cannot describe the entire 
region. This can result in up to ±35% uncertainties in 
column abundances, mostly through uncertainty in the cali- 
bration factor. Calibration for groundbased observations 
assume that the mercurian surface is of homogeneous 
roughness in the visible and near infrared wavelengths and 
can be described by one set of Hapke parameters. We do 
not have the necessary data to know if this assumption is 
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FIG. 4. Disk-resolved brightness curves along the photometric equator (top) and photometric meridian latitude of 50° (bottom) of an orange 
fitter Mariner 10 image. The open diamonds represent the Mariner 10 data. The values of the parameter solution sets are listed in Table 1. 


valid for Na and K wavelengths; it certainly is not true in 
the centimeter wavelength range. 
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